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We have analyzed the spectral density of fluctuations of the energy flux through a mesoscopic 
constriction between two equilibrium reservoirs. It is shown that at finite frequencies, the fluctuating 
energy flux is not related to the thermal conductance of the constriction by the standard fluctuation- 
dissipation theorem, but contains additional noise. The main physical consequence of this extra 
noise is that the fluctuations do not vanish at zero temperature together with the vanishing thermal 
conductance. 
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Fluctuation-dissipation theorem (FDT) [l[ relates the 
fluctuations of a dynamic variable generated at angular 
frequency u> by an equilibrium statistical-mechanical sys- 
tem, to the dissipative part of the response function of 
this system to the force conjugate to this variable at the 
same frequency. This theorem is one of the main physical 
predictions of the "linear response theory" of equilibrium 
transport properties in statistical mechanics, and finds 
applications in practically all areas of the condensed- 
matter physics. The best known case of this theorem is 
the relation between the electric conductance of a resistor 
expressed through Kubo formula Q and the current noise 
generated by this resistor. This relation played important 
role in understanding macroscopic quantum dynamics of 
superconducting structures, where the Josephson effect 
provides a way of directly observing the quantum part 
of current noise as predicted by the FDT [3(. More re- 
cently, this theorem found applications within the studies 
of the mechanisms of decoherence and noise in supercon- 
ducting qubits (see, e.g., Q). From the perspective of 
the linear-response theory, thermal transport represents 
a somewhat special case, since temperature T (more pre- 
cisely, the temperature gradient) that acts as the force 
conjugate to the heat current, does not correspond mi- 
croscopically to any dynamic degree of freedom of a sta- 
tistical system. Nevertheless, it is frequently assumed 
that the FDT holds also for heat transport, and relates 
the spectral density of the energy current S(lo) and the 
heat conductance G t h{^) by (see, e.g., 

S(u) = hujTReG th {uj)coth{huj/2T) , (1) 

This assumption is supported by the fact that at zero 
frequency, Eq. (jTJ) reduces to the form 

2T 2 G th (0), 



S(0) 



(2) 



that follows directly from the fundamental thermody- 
namic result for the magnitude of energy fluctuations. 
Here and everywhere below we define temperature T in 
energy units, by setting = 1. 

The purpose of this work is to demonstrate explicitly 
that the thermal FDT of Eq. Jl) is not valid at non- 
vanishing frequencies, at least as a general statement. 
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FIG. 1: Schematics of a generic constriction admitting energy 
current J between two equilibrium reservoirs with tempera- 
tures Tj and chemical potentials fij, j — 1,2. 



To do this, we calculate the spectral density S(u>) of the 
fluctuations J = J — (J) of the energy current J: 



S(u) 



die 



f [J»J(0) + J(0)J(t)]/2, (3) 



in the set-up characteristic for a "mesoscopic" measure- 
ment of heat transport, see Fig. [T] We consider two 
reservoirs, equilibrated to, in general, different temper- 
atures T\,%, and weakly coupled by a constriction that 
supports propagation of excitations that carry energy be- 
tween the reservoirs. Such a general set-up describes both 
the heat transport by phonons, e.g. in demonstrations of 
the quantization of the phonon heat conductance [6[ , and 
also various structures of "on-chip cryogenics" Ml where 
the heat is transported by electrons or photons [8|,l9|. The 
special status of the temperature T as a parameter in the 
density matrix of the system, instead of being a dynamic 
variable, creates a problem for the microscopic definition 
of thermal conductance. The system Hamiltonian should 
have the property of maintaining local thermal equilib- 
rium required for the temperature and the corresponding 
thermal bias to be well-defined, while at the same time 
creating the non-vanishing energy flux through the sys- 
tem (see, e.g., the discussion in [10j). An important fea- 
ture of the mesoscopic set-up considered here (Fig. 1), is 
that this problem is resolved naturally by clear separation 
of the two processes. Temperature-defining equilibrium 
is maintained inside the reservoirs, while the energy flux 
is created by weak perturbative coupling between them. 
Quantitatively, we first treat the case of phonon heat 
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transport. (With some minor modifications, the same 
considerations and results apply, obviously, to the pho- 
ton heat conduction.) We are interested in the regime of 
frequencies/energies small on the scale set by the phonon 
propagation time along the constriction, u), T/h <C l/r tr , 
where r tr = l/v, with the sound velocity v and constric- 
tion length I. In this low-frequency limit, the short con- 
striction does not lead to any reflection resulting in bal- 
listic propagation of phonons. The simplest description 
of the energy transport by such ballistic phonons in one 
traversal mode starts with the usual Hamiltonian H of 
the field cj>(x) of the longitudinal one-dimensional (ID) 
phonons: 

H = Jdxh(x), h(x) = ^{pl^x)} 2 + k[4>' (x)} 2 }, (4) 

where p and k are the density and compressibility of the 
constriction, so that v — (n/p) 1 / 2 . Writing the Heisen- 
berg equation of motion for the energy density h in the 
form of the continuity equation, h + J' = 0, one obtains 
the following expression for the operator of the energy 
flux J carried by phonons (see, e.g., [Ill]): 



J 



The usual mode expansion of the field <fi, 

^)-(^) 1/2 E(^ eifel - + h-c)/"t /2 , 
iLp k 

where cj^ = v\k\ and L is a normalization length, gives 
J = W E(wfc^p) 1/2 sign(fc)(a fe - 4) (4 ~ °p) ' ( 5 ) 
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(Since none of the quantities we consider below depend 
on x in the small- r tr limit after thermal averaging, we 
set x — in this expression.) Averaging J over the equi- 
librium states of the phonon modes, and taking the ap- 
propriate limit L — > co, one obtains the average energy 
current: 



(j) = 



where rij(E) is the Bose distribution at temperature Tj. 
Expansion of Eq. ([6]) in small temperature difference 
ST/2: Ti,2 = T±5T/2, gives the standard expression for 
the "quantum" of the heat conductance of one phonon 
channel, Gth = TrT/6h. The point of our derivation here 
is that Eq. ([5]) remains valid even if the temperatures 
Ti,2 of the phonons incident on the constriction vary in 
time with frequencies less than l/rt r - This means that 
the resulting expression for the heat conductance Gth is 
also valid for all frequencies in this range. 

Next, substituting the operator J given by Eq. into 
Eq. ([3]) and repeating the same steps that lead to the 



average energy current ((6|), we obtain after some algebra 
the following expression for the spectral density of the 
energy flux noise: 



5 H = ^E/ dEE(E±Q)(l 



n j {E))n j {E±Q) , 



VL = Hlu. This expression generalizes to finite frequencies 
previous results for the noise in phonon heat transport 
(see, e.g., [HI]). Taking the integral, we get 



S(w) 



1 
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[(2^ +n 2 pcothA. (7) 



One can see that even in equilibrium, T\ = T2 = T, the 
spectral density ([7]) is different at fl ^ from the one 
predicted by the FDT, see Eq. ([1]), with the conductance 
G th (u) = TiT/QTi. The FDT reproduces only the first 
part of Eq. ([7]) that corresponds to the T 2 -term in the 
brackets. In addition to this, the full result ([7]) contains 
the fi 2 -term that is no n- vanishing, S{w) = ?i 2 |w| 3 /(247r), 
even at T = 0, when the heat conductance is zero. Phys- 
ically, the origin of this extra term can be traced back 
to finite coupling between the reservoirs, which creates 
quantum fluctuations of their energy even at T = 0, when 
the thermal conductance in the FDT relation (JJJ) vanishes 
since there are no real excitations that could irreversibly 
transfer energy between the reservoirs. In this respect, 
the violation of the FDT for the heat transport discussed 
in this work has the same origin as several other thermo- 
dynamic effects of finite relaxation energy that have been 
discussed in the literature [I3-[l6j]. 

This qualitative picture implies that the breakdown of 
the FDT for heat transport is not a specific feature of 
the phonon heat conduction, but is quite general. To 
demonstrate this, we consider a similar set-up of two 
weakly-coupled reservoirs but in the situation when the 
heat conductance is due to electron propagation between 
them. The main difference with the phonon case is that 
the reflection in the constriction can be non-negligible for 
electrons even when the traversal time r tr is very short on 
the scale set by other energies in the problem. The calcu- 
lation of the finite-frequency heat transport by electrons 
follows the same steps as in the case of phonons. The 
operator of energy density h{x) of ID electrons can be 
written in terms of the electron field ip{ x ) an d the single- 
particle Hamiltonian h = — (K 2 /2m)d 2 / dx 2 + V(x) as 



h(x) = [^hi> + (h^)ip}/2. 



(8) 



The symmetrized expression ([8]) is needed to ensure that 
h(x) is Hermitian. Then, the Heisenberg equation of mo- 
tion for h{x) with the Hamiltonian H = J dxh(x) takes 
the form of the continuity equation, h + J' = 0, with the 
energy flux operator 



J = {-in/Am)\^{h^)' - {^)'h^ - h.c] 



(9) 
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This expression shows that if one decomposes the field ip 
into the stationary scattering modes, the energy current 
has the same form as the usual probability current, the 
only difference being that h multiplies each mode by its 
energy. Explicitly, introducing the creation/annihilation 
amplitudes a\,a,k for electrons incident from one elec- 
trode, and b\, bk from the other, we obtain the following 
mode expansion of J: 

J= V fJ2 e -^[D(ala p ~blb p ) + Vmi(a{b p + bla p )} . 

(10) 

Here vf is the Fermi velocity, D and R are, respectively, 
the transmission and reflection probabilities of the con- 
striction, D + R = 1, and ek lP are the electron ener- 
gies. Also, we assume that both the Fermi energies in 
the electrodes, and h/T tr are much larger than the typi- 
cal excitation energies T, SI, and eV, where V is the bias 
voltage between the electrodes. This implies, in partic- 
ular, that the scattering probabilities D, R are constant 
in the energy range of interest. The energy current J 
(fTOf corresponds directly to the heat flow into/out of 
the reservoir j, if the electron energies ek, P are measured 
in Eq. (fTOf relative to the chemical potential of this 
reservoir. (Note that in the phonon calculation above, 
this condition was satisfied automatically, since /i = 
for phonons.) If V = 0, so that there is no shift between 
the chemical potentials of the two electrodes, the aver- 
age of Eq. (fT0|) represents both the heat flow, —J\, out of 
one electrode and the heat flow, J2, into the other one: 
— (Ji) = (J2) = (J)- The total average generated heat is 
zero, (J1 + J2) — 0. If, however, V 7^ 0, then one needs to 
measure the energies relative to the two different levels [i\ 
and /i2 in the two reservoirs, and — (Ji) 7^ (Ji)- The dif- 
ference between the two heat flows is obtained by replac- 
ing the energy (e^ + e p )/2 in Eq. (|10p with ^1—^2 = e V- 
After this substitution, Eq. (fTOf reduces to IV, where I 
is the electric current between the reservoirs. Therefore, 
in the case of non- vanishing bias voltage V, the total gen- 
erated heat is non- vanishing, and equal to the Joule heat, 
( Ji + J2) = (I) V, where the individual heat flows Jj are 
obtained from Eq. (|10p by measuring the energies ek, P rel- 
ative to the chemical potential fij of the corresponding 
electrode. 

Taking thermal average and the limit L — > 00 in 
Eq. (pTO]) . we find the average heat currents (Jj), j = 1,2, 
into the two electrodes: 

(Jj) = / dE(E - fij)[fi(E) - f 2 (E)] , 

where fj(E) is the Fermi distribution function of elec- 
trons in the jth reservoir. This gives for the heat cur- 
rents: 



The two terms in this expression represent, respectively, 
the usual Joule heating, which in this case is distributed 
equally between the electrodes, and the heat transport 
between them. For small temperature difference between 
the electrodes, Eq. (fTTj) gives the thermal conductance, 
Gth — kDT/QTi, that coincides with the phonon conduc- 
tance. As with the phonons, an important point of our 
derivation here is that this thermal conductance is inde- 
pendent of frequency in the considered frequency range 
below 1 /r tr and the frequencies set by the Fermi energies 
in the electrodes, and the energy scale of the variations 
of the transmission probability. 

At V 7^ 0, when the total heat fluxes into the two 
electrodes are different due to Joule heating, to describe 
specifically the heat transfer between the electrodes, one 
needs to define the heat current as J = (J2 — Ji)/2. 
Indeed, as one can see from Eq. (1TT1) . the average heat 
current (J) defined this way is not affected by the Joule 
heating. This definition corresponds to the simple pre- 
scription of measuring all energies in Eq. (|10|) relative to 
the midpoint between the chemical potentials of the two 
electrodes. Then, the same steps as for the average cur- 
rent, including thermal averaging and the L — > 00 limit, 
give for the spectral density of the heat current noise 

S ^ = S / dE[E ± 0/2)2 ^ {E) ' 

(1 - fj(E ± n)) + Rfj(E)(l + f r (E ± n))] , (12) 

where j' is defined as j' 7^ j with j,j' = 1,2. One can 
see that without the bias voltage and reflection, V = 0, 
R = 0, Eq. (fT2|) coincides with the phonon result (0. For 
equal temperatures of the electrodes, T\ — T2, calcula- 
tion of the integral (fT2|) gives 

S(oj) = -J^—{2D\n 2 + (2ttT) 2 +3(eV0 2 lftcoth — 

+Rj2[(eV ± nf + (2irT) 2 ](eV ± SI) coth eV Jl - }.(13) 

This result (as well as Eq. ([7]) for the phonons) can be 
extended naturally to the situation when the contact be- 
tween the electrodes supports many electron modes with 
transparencies Dk- In the equilibrium case that is of the 
main interest here, Eq. (|13l) gives the following expression 
for the heat noise in such a multi-mode contact: 

S(lu) = (G/12e 2 )[fi 2 + (2irT) 2 ]Slcoth(Sl/2T) . (14) 

Here G = J2k Dke 2 / (2irh) is the electric conductance of 
the contact, which is related by the Wiedemann- Franz 
law to the heat conductance G t h = tt 2 GT '/ '(3e 2 ). 

One of the interesting physics features of Eq. ([T3|) in 
the case of non-vanishing V and R is the shot noise of 
the energy current associated with the individual elec- 
tron scattering events. Energy currents carried by indi- 
vidual electrons were predicted in [l7| and demonstrated 
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recently in the form of RF-cooling in a metallic single- 
electron transistor with an alternating voltage at the gate 
fl8j |. Equation (fT3l shows that in a biased contact, the 
average heat current is accompanied by shot noise of heat 
due to the scattering of discrete electrons. Quantita- 
tively, the same multi-mode generalization of (|13p gives 
this noise for \eV\ T,fl as 

S(V) = eFG\V\ 3 /12, (15) 

where F = J^k ~ Dk)/J2 k is the standard Fano 
factor that characterizes the shot noise of electric current. 

Returning to Eq. (fT4")) and quantum fluctuations of 
heat, we see again that, similarly to the situation with the 
phonon heat conductance, equilibrium fluctuations of the 
heat current are not described correctly by FDT at finite 
frequencies. While the heat conductance of the contact 
is frequency- independent in the range discussed above, 
the noise contains frequency-dependent part which does 
not vanish at T = together with the heat conduc- 
tance. Physically, these fluctuations are produced by 
virtual electron transitions between the two electrodes 
due to finite coupling between them. This mechanism 
is the same as for the quantum fluctuations of electrical 
current in the contact, and Eq. (fl"3f for the energy fluc- 
tuations is quite similar to the corresponding equation 
for the current fluctuations. The fact that the quantum 
fluctuations of electric current are still consistent with 
the FDT, in particular, vanish together with the electric 
conductance G, while the thermal fluctuations do not 
agree with FDT, is a reflection of the special nature of 
temperature in statistical mechanics. 

As the last point of our discussion of the thermal FDT, 
we would like to make more explicit the set of assump- 
tions underlying the notion of the frequency-dependent 
thermal conductance G t h{u). We do this in the case 
of electron heat transport considered above, limiting the 
discussion to the tunnel approximation D <C 1, when 
electron scattering can be described with the usual tun- 
nel Hamiltonian 

H = H 1 +H 2 + H T , H T =J2( T kp4 b P + h - c -)i ( 16 ) 

k,p 

which explicitly separates Hamiltonian of the electrodes 
H\ 2 and the tunneling term Ht- Here Tk p are the tunnel- 
ing amplitudes that can be expressed through the junc- 
tion conductance G. The heat flow defining the thermal 
conductance at finite frequencies should be driven by a 
small time-dependent temperature difference 5T(t). Sep- 
arating one frequency oj, we take Ti j2 = T± (ST/2)e~ tult . 
Expansion of the density matrix of the equilibrium elec- 
trodes in ST <C T gives the <5T-induced correction Sp to 
it as 

Sp{t) = -p Q Q{5T/2T 2 )e- wt . (17) 

Here po = (l/Z)e~( Hl+H2 " T is the main part of the den- 
sity matrix and Q = Q' — (Q'), with Q' = H 2 — Hi 



and (...) denoting the average over po- We note that 
the assumption of equilibrium with the time-dependent 
temperature requires that lo <C 1/te, where te is the 
relaxation time in the electrodes. This condition is not 
satisfied directly in the model of non-interacting phonons 
or electrons for which we calculated the energy flux noise 
S(ui). One can still use Eq. (IT7l) for the non-interacting 
particles incident on the junction, assuming that the 
temperature-defining relaxation is concentrated infinitely 
deep inside the electrodes. The usual perturbation the- 
ory in Ht around Sp (|17p then gives 

G th (w) = ^ J™ dte iut (J(t)H T Q - QH T J(t)), (18) 
where the energy flux operator is 

J = Q/2 = (i/2ft) 5^(c fc + e p ){T kp a\b p - h.c.) . 

k,p 

Equation for the thermal conductance similar to (fT8|) 
can also be derived in a general situation, when the heat 
current is driven by some arbitrary relaxation interac- 
tion V (and not the tunneling Ht)- Formally, the break- 
down of the FDT for the heat transport discussed in this 
work arises from the difference between the structure 
of Eq. (fl8|) and expressions for the "dynamic" linear- 
response coefficients, e.g., the electric conductance. Ex- 
plicitly, evaluating JT8]) we get ReG th (u) = 7r 2 GT/(3e 2 ), 
which is frequency-independent in agreement with the 
arguments provided above, and in contradiction to the 
FDT (JlJ , if compared with the energy flux noise ([M)) . 
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